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Abstract 

We calculate the masses of the low-lying states with quantum numbers J PC = 
++ , 1 in the Higgs and confinement regions of the three-dimensional SU(2) 
Higgs model, which plays an important role in the description of the thermody- 
namic properties of the standard model at finite temperatures. We extract the 
masses from correlation functions of gauge-invariant operators which are calcu- 
lated by means of a lattice Monte Carlo simulation. The projection properties 
of our lattice operators onto the lowest states are greatly improved by the use of 
smearing techniques. We also consider cross correlations between various operators 
with the same quantum numbers. Prom these the mass eigenstates are determined 
by means of a variational calculation. In the symmetric phase, we find that some of 
the ground state masses are about 30% lighter than those reported from previous 
simulations. We also obtain the masses of the first few excited states in the sym- 
metric phase. Remarkable among these is the occurrence of a ++ state composed 
almost entirely of gauge degrees of freedom. The mass of this state, as well as that 
of its first excitations, is nearly identical to the corresponding glueball states in 
three-dimensional SU(2) pure gauge theory, indicating an approximate decoupling 
of the pure gauge sector from the Higgs sector of the model. We perform a de- 
tailed study of finite size effects and extrapolate the lattice mass spectrum to the 
continuum. 



1 Introduction 



The study of three-dimensional field theories has attracted a lot of attention over the 
past few years. While some models are investigated for field theoretic reasons or because 
they are more easily accessible than their four-dimensional homologues, others have an 
immediate physical meaning in the context of four- dimensional field theory at finite 
temperature. It has been known for a long time that for temperatures much higher than 
any mass scale of a given theory its non-static Matsubara modes may be integrated out 
perturbatively to yield a three-dimensional effective theory for the zero modes This 
effective theory describes the static long-range physics of the underlying four-dimensional 
finite-temperature theory, and moreover contains all the infrared divergences and non- 
perturbative phenomena that spoil a purely perturbative treatment of the latter. 

In particular, the three-dimensional SU(2) Higgs model represents an effective high- 
temperature theory for the standard electroweak model, after neglecting the U(l) sector 
and fermions in a first approximation. Since it was realised that the baryon asymmetry of 
the universe could have been generated during a first-order electroweak phase transition 
[0, a lot of effort has been spent to determine the order and dynamics of this transition 
in detail. The perturbative procedure of dimensionally reducing the four-dimensional 
SU(2) Higgs model at finite temperature to a three-dimensional effective theory has 
been carried out in great detail in ref. |§, and the corresponding relations between the 
three-dimensional and four- dimensional parameters and temperature may be found there. 
Bearing these relations in mind, we shall stay entirely within the framework of the 2+1 
dimensional SU(2) Higgs model in this paper. 

There are already several analytical Q-0 and numerical |§-|l(| studies of the three- 
dimensional SU(2) Higgs model. While the main motivation for these studies was the 
phase transition itself, there also emerged the problem of understanding the structure 
of the symmetric phase. Due to infrared divergences in vector loops, straightforward 
perturbation theory breaks down in the symmetric phase, and until recently not much 
was known about the particle spectrum and the effective interactions in this parameter 
region of the theory. In a recent analytic calculation, the loop expansion was reorganised 
by resumming masses and vertices, which led to a set of gauge independent gap equa- 
tions for the vector boson and the Higgs masses, defined on their respective mass shells. 
The solutions of these equations predict a non-vanishing vector boson mass and scalar 
vacuum expectation value in the "symmetric phase" , which thus would be interpreted as 
another Higgs phase, just with different parameters. On the other hand, lattice simula- 
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tions in four dimensions at finite temperature ['1 1 as well as in three dimensions ||, 10 1 
have reported vector boson masses about four times larger than predicted by the gap 
equations.^] The picture conjectured from the lattice results is one of a symmetric phase 
with confining behaviour (in the sense of QCD) and a particle spectrum consisting of 
bound states. Similar conclusions may be drawn from analytic investigations of truncated 
renormalisation group equations which indicate strong coupling effects in the symmet- 
ric phase 0. The picture of a QCD-like symmetric electroweak phase has also been 
employed for a model calculation of bound state masses ||. One possible explanation 
of the large discrepancy between the two approaches is a breakdown of the resummed 
loop expansion of H in higher orders. In principle, however, it is also conceivable that 
an extremely low-lying state might not have been visible on the lattice sizes that have 
been investigated. Moreover, all simulations have emphasised the difficulty of measuring 
correlation functions in the symmetric phase due to the extremely low signal-to-noise 
ratio. 

The purpose of the present paper is to shed more light on the situation in the sym- 
metric phase by employing new techniques which allow a more reliable lattice calculation 
of the mass spectrum. The masses are extracted from correlation functions of gauge in- 
variant operators. In order to improve the sensitivity to low-lying or bound states we 
construct a large set of non-local operators by employing a "blocking" technique simi- 
lar to the one which has proved to be useful in pure gauge calculations |L2], [Tj|. As we 
shall see, this procedure turns out to be very effective in enhancing the projection of our 
operators onto the lowest states. Moreover, it reduces the statistical errors significantly, 
yielding more accurate results for the masses. We also measure cross correlations be- 
tween different operators. Diagonalisation of the corresponding correlation matrix then 
unmixes the superposition of the ground state and the excited states. This procedure 
further improves the signal for the lowest states and, more importantly, enables us to 
estimate the masses of the first few excited states. It also allows us to extract informa- 
tion about the overlap any individual operator has with a given state, and hence the 
coupling between the different states. We then perform a study of finite-size effects and 
an extrapolation to the continuum limit for two points in parameter space, one located 
in the symmetric and one in the Higgs phase. 

The paper is organised as follows. In section |2] the lattice action and the basic oper- 
ators used in mass calculations are discussed. The details and more technical aspects of 



1 For a detailed comparison of the lattice and analytic approaches in three and four dimensions in 
the context of the electroweak phase transition, see || . 
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our simulation are described in section |3|. In section^ we present our results, analysing 
in detail the reliability of our mass estimates of the ground state, and including the 
extrapolation to the continuum limit. Finally, section^ contains our conclusions. 



2 Action and basic operators 

The action of the SU(2) Higgs model in 2+1 dimensions and its general properties in the 
continuum and on the lattice have been discussed previously in the literature f3H-||10||. 
In order to fix the notation and to give all equations used in this paper we list some of 
these general aspects here. 

The continuum action is given by 



S 



d 3 x Tr 



where all fields are in a 2 x 2 matrix notation, 



(1) 
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The gauge coupling g 3 and the scalar coupling A3 have mass dimension 1/2 and 1, 
respectively. The action ([I]) can be parametrised by two dimensionless parameters, which 
may be chosen to be X^/g^ and fJ^/ 93- Fixing these parameters determines the physical 
properties of the theory. The corresponding lattice action may be defined as 

S[U,4>] = /3 G E(l-^Trf/ p )+E{-^E^Tr(0t( x )^( x )0( x + ; a 

p V Z J x { u=l 1 



+^Tr( ( j ) \x) ( j ) (x)) +(3 R [iTr(0t(x)0(a;)) - lj" 



(3) 



Due to the super-renormalisability of the theory (JJ), the only parameter receiving 
ultraviolet renormalisation is the scalar mass parameter n^/g^, whose corrections have 
been determined at the two-loop level in perturbation theory using the MS scheme. 
The corresponding two-loop calculation in lattice perturbation theory was carried out in 
nj . Requiring that the renormalised mass parameters be the same in both regular isat ion 



schemes, a relation between the parameters labelling the continuum and lattice theories 
has been established ilHI, 
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(6) 



16 ' ~9l \9l) ) \ 2 V ■ g$ 

with the numerical constants X = 3.17591 and ( = 0.09. 

A Monte Carlo simulation of any quantity in the theory (0) is carried out for a 
given set of bare parameters Pq,Ph,Pr- In order to establish contact with the desired 
continuum physics one first has to perform an infinite volume limit, i.e., do simulations 
on lattices much larger than the largest correlation length of the theory such that the 
results do not show any dependence on the lattice size. Secondly, one has to perform 
a continuum limit a — > 0, i.e., to simulate at different values of (3q and extrapolate to 
(3g — > oo. This limit has to be taken in such a way that the renormalised quantities 
parametrizing the theory remain constant. The corresponding 'lines of constant physics' 
in the space of the lattice parameters are given by equations @-(U). 

The actions (|l|) and (H) have an SU(2)i ocal x SU(2) globa i symmetry. Physical states 
are described by gauge invariant operators. After decomposing <p(x) as 

<P(x) =p(x)a{x), p 2 (x) = ^Tr , p(x) > 0, a(x) E SU (2) , (7) 

one may define the gauge-invariant composite field 

V^( x ) = ar (x)Un(x)ac(x + p) . (8) 

While p(x) and V^x) are invariant under local transformations, they transform under 
the diagonal global SU(2)di ag subgroup, customarily termed weak isospin, as 

p'(x) = p(x), V^x) = AV^x)A-\ AGSU(2) diag , (9) 

i.e., the lowest excitation of p(x) describes the isoscalar Higgs boson while the matrix- 
valued Vp(x) transforms as an isovector. A single field representing the spin-one, isospin- 
one W boson may be obtained from the composite link variable by taking the trace with 
an insertion of a Pauli matrix r a . Taking the trace without reinsertion produces another 
spin zero isoscalar operator. A third ++ isoscalar operator is given by the plaquette. 
Thus we consider the following set of basic operators for the description of the low-lying 
states, 

++ : R(x) = ^Tr ^{x)<f){x] 
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++ : L(x) = J2 o Tr (a\x)U^(x)a(x + /})) , 

++ : P(x) = U 1 (x)U 2 (x + i)Ul(x + 2)U 1 2 (x), 
1— : V;(x) = ^Tv (r a a\x)U^x)a(x + fl)) . (10) 

The plaquette P is particularly interesting because it consists of gauge degrees of freedom 
only, and in the three-dimensional pure gauge theory it is the simplest operator one can 
use to describe the ++ glueball In the theory with scalars one expects a mixing 
of this operator with the other ++ operators due to the coupling of gauge and scalar 
degrees of freedom. 

The phase structure of the three-dimensional model with lattice action (|3D has not 
been fully mapped out by numerical simulations as yet. However, from analyticity consid- 
erations |T3| and numerical studies of the four-dimensional model |]16| one expects the fol- 



lowing qualitative picture: the three-dimensional parameter space spanned by j3c, Ph, Pr 
is divided into Higgs and confinement-like (or symmetric) regions by a surface of first- 
order phase transitions which is crossed by changing (3h for fixed /3r- At sufficiently 
large values of (3r and small values of fie this surface is expected to terminate so that 
the two regions are analytically connected. In this region, there is no phase transition 
but just a crossover as (3h is varied. Numerically, however, this region in the phase dia- 
gram has so far not been accessed in the three-dimensional theory. The continuum limit 
is represented by a single point in the phase diagram, (3g — * oo, (3r — > 0,(3h — ► 1/3. In 
order to describe different continuum theories, the continuum limit has to be taken along 
different paths in the parameter space, as specified by equations (f|)-(|])- 

The term confinement-like is chosen to distinguish the behaviour of the theory in this 
region of parameter space from the confinement realised in the three-dimensional pure 
gauge theory. There the potential between static charges rises linearly with distance, 
without any bound. In the Higgs model one expects a flattening of the potential at some 
large distance, due to pair creation of scalars breaking the string between the static 
charges, just as fermions break the string in QCD. From the analytic connectedness of 
the Higgs and the confinement regions it follows that for every state in the Higgs region, 
there is a corresponding one in the confinement region. In particular, the same operators 
fllPf ) may be used to describe physical states in both regimes. The global isospin symmetry 
is realised in the Higgs as well in the confinement region, so one expects low-lying states 
with the same quantum numbers in both regions. 
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3 The Simulation 



The purpose of this paper is a closer investigation of the mass spectrum in the 
confinement-like phase, by calculating correlations of operators of the type in eq. fllPp. 

In this section, we describe the details of our calculation, including the simulation 
algorithm, the construction of "blocked" or "smeared" operators for the Higgs and vector 
bosons and the way in which these are used to compute matrix correlators. We conclude 
this section with details about the statistical analysis and the fitting procedure employed 
to obtain our final mass estimates. 



3.1 Simulation algorithm and parameters 

Our Monte Carlo simulation is performed using the lattice action in eq. (|3]), containing 
the bare parameters /3q, (3h and (3r. 

For the update of the gauge variables we use a combination of the standard heatbath 
and over- relaxation algorithms for SU(2) [|T7| , 18 1. The scalar degrees of freedom are up- 



dated using the algorithm proposed in |Tj|, which uses the four real components of the 
scalar field 4>(x). Thus, no separate updates of the radial and angular parts p(x),a(x) 
are required, leading to a simple implementation of the algorithm. As explained in [|11J , 
over-relaxation (reflection) steps in the update of the scalar field can be easily incorpo- 
rated, provided the Higgs self-coupling (3r is not too large, which would lead to a poor 
acceptance rate. In our simulation, where {3r = O(10~ 4 ), we achieved acceptance rates 
of well over 90%. Higher values of [3r could be simulated, for instance, by using the 
reflection algorithm described in ||20|| . 

In our simulation, a "compound" sweep consists of a combination of heatbath (HB) 
and several reflection (REF) updates of the gauge and scalar fields, 

1 HB{[/} + 1 HB{0} + n OR [REF{U} + n rcf REF{0}}. (11) 



In accordance with ref . |19| , we chose uor to be roughly equal to the inverse scalar mass 
in order to achieve maximum decorrelation. With this choice we found that the average 
integrated autocorrelation time estimated using the scalar mass was close to one, in 



agreement with [plj 



Our simulations were performed for inverse gauge couplings (3g = 7, 9 and 12. We 
restricted our attention to one point in the symmetric and one point in the Higgs region 
of parameter space chosen sufficiently away from the phase transition, so that the system 
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does not tunnel between the phases. In order to compare our results directly with those 
of a previous calculation of the lightest scalar and vector masses 0, we work at the same 
fixed value of A 3/^3, 

A3 = /3r(3g = Q (12) 
9s Ph 

which in the context of the four-dimensional theory corresponds to a Higgs mass at 
tree-level and zero temperature of Mh — 35 GeV. 

In the symmetric phase of the model we initially chose j3c = 12, (3h = 0.3411. The 
value of (3 n was then fixed by the relations (|J)-(§|). The same relations determine the 
continuum scalar mass parameter in the MS scheme as n\/g% = 0.089. Our point in 
the Higgs phase of the model was fixed to be (5q = 12, (3h = 0.3418, which corresponds 
to i&,jg\ = —0.020 in the continuum. At (5q — 7 and 9 the corresponding values of (3h 
and (3r were chosen according to the "lines of constant physics", eqs. (§])-©, using the 
constraint eq. fli~2|) . 

At /3q = 12, Monte Carlo runs were performed on lattice sizes ranging from 10 2 • 12 
up to 40 3 in order to analyse finite-size effects in detail. This is of special importance in 
the symmetric phase of the model, where we are particularly interested in the possible 
occurence of very light states. 

For all our observables, statistics were gathered from about 30 000 compound sweeps. 
In a few cases, statistics were increased to a total of 75 000 sweeps. 



3.2 Constructing improved operators 

The main difficulty encountered in recent attempts to compute the mass spectrum in 
the symmetric phase of the SU(2) Higgs model JO], |J, was the low signal-to-noise ratio 
in the computation of the correlation function 

C(t) = Ee ip - (x - x V(x,tMx',0)> c = £ \(0\cp(0)\n)\ 2 e- aEnt (13) 

x,x' n>0 

= £ |c n | 2 e- aB ^~~| Cl | 2 e- aBl <, (14) 

n>0 

where <£>(x, t) denotes any one of the operators in eqs. fllOD, and E n > i£ n _i is implied. For 
our numerical calculation of the masses we use the zero momentum timesclice averages 
of the original operators, i.e., p = in the above expression. 

One important goal of our simulations is to investigate the possible existence of 
very low-lying states in the symmetric phase of the model, such as predicted by the 
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analytic approach in ||]. If an operator has a bad projection onto the lightest state 
one must be able to follow the signal to sufficiently large t before the ground state 
dominates C(t). A poor signal-to- noise ratio of the correlation function will then hamper 
any effort to establish the existence of such a state. The problem is further exacerbated 
in the symmetric phase, where, due to the confining behaviour of the theory, the particle 
spectrum may consist of bound states, having a larger spatial extension than their point- 
like counterparts in the Higgs region. Previous experience with calculations of the glueball 
spectrum in pure gauge theory shows that conventional local operators have indeed a 



bad projection onto bound states in confining theories [12], |13| . The situation could 
be considerably improved by constructing "blocked" or "smeared", non-local operators 
[fL2j which are of similar extended structure as the bound states they are supposed 
to project on. Similar techniques, which preserve gauge invariance, have been developed 
and successfully applied in simulations of lattice QCD [p2|. 



Here we are applying and extending these ideas in order to construct non-local ver- 
sions of the operators defined in eq. fllPP. Some of these techniques were applied to the 
four-dimensional SU(2) Higgs model in ref. p3 . 



Following [|T^| , we construct composite ( "blocked" ) link variables Ujfi (x) of blocking 



level n according to 

U { ;\x) = U^ l \x)U^ l \x + fi) (15) 

±2 

+ E U^ n - l \x)U^- l \x + v)U^- x \x + + v)U^ n - l) \x + 2 ft) . 

The links U ( jj l \x) are twice as long as those at the lower blocking level n — 1. We shall 
refer to this as "link blocking" in the following. It seems natural to design a similar 
procedure for the scalar fields 4>(x). A "site-blocked" scalar field (j)^(x) at blocking level 
n can be constructed iteratively from a field at a given lattice site and its covariant 
connection with the four nearest neighbours, 

2 

0W( X ) = ^-^(x) + [U^ix^^ix + p,) + UjT^Hx - £)0 (n_1) (z - fi)] ■ (16) 

Ai=l 

Clearly both blocking procedures can be iterated, thereby quickly increasing the number 
of links and sites contributing to a given composite variable. 

Non-local blocked operators are now constructed from the basic ones ([H]) by replac- 
ing the scalar and link variables with composite ones at a desired blocking level. Note 
that the blocking steps are constructed in a way which preserves the gauge invariance of 
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the original operators. The basic operators R and P contain only site and link variables, 
respectively. By applying the corresponding blocking procedure to these operators, we 
get N operators of different spatial extension, where N denotes the maximal blocking 
level. These we write as PS n \x) and P^ n \ with n = 0,...N. For the operators L and 
V both site and link blocking can be applied, so from each of them we construct a set 
of N x N operators, denoted by L^ m \x) and V^ nm \x), n, m = 0, ...N, where the first 
upper index stands for site and the second for link blocking. 

3.3 Cross correlations 

The blocking procedure described in the previous subsection is designed to yield an 
optimal operator for a given set of quantum numbers. In an attempt to further separate 
the excitations from the ground state we can also utilise the information contained in 
our non-optimal operators by considering cross correlations between different operators 
in the same channel. 

For a given set of quantum numbers we construct a set of, say, N lattice operators, 
fa : i = 1,..,N, with those quantum numbers. We normalise these operators so that 
(fa fa) = 1, and we impose the same normalisation on all the operators we discuss 
below. To find the energy of the lightest state we use a variational criterion. That is to 
say, we find the linear combination of the fa that maximises 

<0t(a)0(O)> = (^e-^) . (17) 

Call this operator $!. In the limit where the basis {fa} is complete, this procedure 
becomes exact. That is to say, if the lightest state is |1) and the corresponding energy is 
Ei, then 

$i|mc) = |1) , (18) 

and 

e -°Ei = ($^(^(0)) . (19) 

We can find higher excited states by a simple extension of this procedure. Let the first 
excitation be \2) and let the corresponding energy be E 2 . We consider the subspace {</>«}' 
of {fa} that is orthogonal to $i, i.e., such that (<Ey (0)0(0)) = 0. We apply the same 
variational criterion as above, but restricted to this subspace. This gives us an operator 
$2- in the limit where our original basis becomes complete, we have 

<& 2 \vac) = |2) , (20) 
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and 

e -aE 2 = ($ 2 t( a )$ 2 ( )) . (21) 

We can continue this procedure obtaining operators <3> 3 , $ 4 , ... from which we can obtain 
the energies of higher excited states. 

In our case our basis is finite, and we can obtain at most N operators With 
such a limited basis eq. ([19]) provides at best an estimate for aE]_. We improve upon this 
estimate by calculating the correlation function ($ 1 ^(/f:)$ 1 (0)) for all t. If we define an 
effective energy by 

e -aE cS (t)t = ($^(^(0)) , (22) 

then we know that E e ^{t) will approach Ei from above as t increases. The more effective 
our variational procedure, the smaller the value of t at which this occurs (for a basis that 
is complete we would find E e g(t) = E\ for all t). So we can estimate aE\ from the value of 
aE c g(t) on its 'plateau'. In practice, what we actually do is to fit the correlation function 
to an exponential in t for t large enough (as described below). From the exponent we 
then obtain our estimate for aE\. From the coefficient of the exponential we obtain the 
normalised projection of our operator onto the lightest state, i.e. |(l|$i|t>ac)| 2 . If we 
have a good basis of operators then this projection will be close to one. In practice this 
is always the case in the scalar channel, where the projection is often consistent with 
unity. In the vector channel the projection tends to be ~ 0.8. 

We follow the same procedure for excited states, extracting aEi by fitting an expo- 
nential to ($j(t)<l ) j(0)) for large enough t. One must be more careful here than with the 
ground state because, with a finite basis, the operator $j will have some projection onto 
all states, not just onto So as t — > oo its correlation function will ultimately vary as 
~ exp(— Ext) and not as ~ exp(— Eit). Thus, by fitting an exponential at larger t we may 
underestimate the value of aEi. In practice this is not a problem where the operators are 
very good. For example if the projection of $i onto the lightest state is 1 — e, then the 
projection of $2 onto this lightest state is < e. If e is as small as it is in our calculations, 
then this potential contamination of E 2 by E\ is insignificant. The same argument can be 
used for higher excited states. In general, where we quote a mass without qualifications, 
we are confident, by examining the relevant projections, that our mass estimate is not 
significantly contaminated by admixtures of any of the lighter states that we list. 

In practice our lattice is finite and so in the above we replace e~ Et by e~ Et + e~ E ^ T ^ 
where T is the length of the lattice in the t-direction. 

The procedure we follow to obtain the $j is standard Define the N x N corre- 



10 



lation matrix C(t) by 

Cij(t) = (^m(0)) . (23) 
Let the eigenvectors of the matrix C _1 (0)C(a) be v l ; i = 1, . . . , N. Then 

N N 

<$>i = CiJ2 v l0k = J2 a ^k , (24) 

k=l k=l 

where the constant Cj is chosen so that is normalised to unity. 

We would like to emphasise that there are many possible variations on the above 
variational procedure. For example we could apply it to t = 2a rather than to t — a. As 
a check we have performed such an alternative calculation. We further remark that, in 
practice, the best of our original <pi operators is already so good that the calculation of 
the ground state in each channel is not greatly improved by going to the <3?j operators. 
It is if we wish to obtain the excited states that this analysis becomes indispensable. 

In our actual calculations, we typically compute a 9 x 9 matrix of correlators in 
the ++ channel, which consists of the three operators R, P and L, each taken at three 
different blocking levels. In the 1 channel, where only operators of type V are known, 
three different blocking levels are used to compute a 3 x 3 correlation matrix. 



3.4 Fits and error analysis 

All our mass estimates are obtained from measured correlation functions of operators 
in the diagonalised basis defined in the preceding subsection. The ansatz we use for the 
asymptotic behaviour of the correlation function on a finite lattice, for large T, is 

Ci(t) = ($}(*)$i(0)) = A (e~ aMit + e- aM ^ T "')) , (25) 

where % labels the operator, and T denotes the extent of the lattice in the time direction. 
This expression would be exact for all t if the basis of operators was complete. To monitor 
deviations from this behaviour we define an effective mass according to 

aM eS (t) = arcoshj ^ + ^+^-^ j , (26) 

where C(t) denotes either Ci(t) or Ca(t). As one readily sees, this definition has the 
desired property that aM eS (t) = aMi for those t where C(t) is accurately given by (|25|). 

Estimates for the masses aMi and amplitudes A, are obtained from correlated fits of 
Ci(t) to eq. (p5|) over a finite interval [ti, t 2 }. Our choice of the fitting interval is guided by 
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the plateaux observed in the effective masses (0), and is constrained by the requirement 
that a reasonable x 2 /dof should be obtained. 

Our individual measurements of Ci(t) are accumulated in bins of typically 500 mea- 
surements each. Statistical errors on aMi and A{ are obtained from a jackknife analysis 
of the fits to the average of Ci(t) in each jackknifed bin. 

It has been known for some time that correlated fits may amplify hidden systematic 



errors in the data p5| . Therefore we repeated all our fits using an uncorrelated covariance 
matrix. The difference between the results obtained using either correlated or uncorre- 
lated fits are quoted as a (symmetric) systematic error on our mass estimates. In most 
cases we found the systematic error arising from this procedure much smaller than the 
statistical error. For the final extrapolation of masses and mass ratios to the continuum 
limit, statistical and systematic errors are added in quadrature before the extrapolation 
is performed. 



4 Results 



In this section, we present our main results. We start with a discussion of the effects of 
the blocking procedure and the diagonalisation of operators in subsections [O and pT2", 
using our data at /3q = 12 on the largest lattices we investigated in the confinement 
phase (40 3 ), and in the Higgs phase (20 3 ). The main results on the spectrum, which 
were obtained using diagonalised operators at all three values of 0q, are presented in 



subsections [4. 3| and |4.4| . Finally, in subsection |4.6| we give our mass estimates extrapolated 
to the continuum limit. 



4.1 The effects of the blocking procedure 

A priori nothing is known about the projection properties of the individual operators 
in our set {<pi}. The candidates with the best projection onto the lowest states have 
to be determined from actual simulations. A criterion to judge the performance of an 
operator is its effective mass at time separation one, where the lowest value indicates 
the least contamination from excited states. Figure p] illustrates the effect of the blocking 
procedure for the purely scalar/gauge operators R/P (cf. eqs.flTOD) in the Higgs phase and 
the confinement phase, respectively. In the Higgs phase, nothing is gained by blocking 
the R operator, while in the confinement phase four iterations are necessary before it 
reaches its optimal projection. For the plaquette P, three blocking steps are required 
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Figure 1: Effect of blocking on the R and P operators in the Higgs (left) and confinement 
(right) phases. 

to get to the minimal effective mass in either phase, but the improvement is far more 
pronounced in the confinement phase. 

Figure ^| shows effective mass plots for the operators with the best projection of each 
basic type in the ++ channel, again for the Higgs and confinement phases, respectively. 
In both regimes the operator with the best projection onto the ground state is of the 
type R, with nearly 100% overlap at the optimal blocking level. The ground state could 
in principle also be extracted from the correlation function of the best candidate of the 
type L at large time separations. However, its projection is much worse, five to six lattice 
spacings are needed until excitations have died away, and a mass calculated from this 
correlation function would be much less accurate. Of particular interest is the behaviour 
of the plaquette correlations. While they are dominated by noise in the Higgs phase, they 
suggest a separate plateau in the effective mass plot in the confinement phase. Up to 
those time separations for which we have a good signal, there seems to be no tendency 
for this operator to mix with the other ++ operators. We shall return to this observation 
below. 

Finally, Figure ^| shows the result of the blocking procedure on the effective masses 
of the vector boson. In the Higgs phase blocking slightly improves the projection of the 
operator V, but it is not difficult to extract a mass also from the unblocked one. In the 
confinement phase the situation is rather different. The unblocked operator does not give 
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Figure 2: T/ie optimally blocked operator of every basic type in the ++ channel, Higgs 
(left) and confinement (right) phase. 

any signal beyond noise, and the figure displays how even for the best blocked candidates 
excitations die out only very slowly. 

In particular, the figure illustrates how one might easily extract too large a mass for 
the vector boson, if one only used a non-optimal operator such as the one symbolised 
by the triangles. The effective masses produced by the different operators do not seem 
to merge at a common ground state up to the distances to which we can follow the 
signal. An indication that we do really see the ground state is the fact that further 
iterations in the blocking procedures for either links or scalars, beyond the level of the 
operator V 55 , again result in a worse projection. Figure ^ also displays a nice side effect 
of the blocking procedure. Since the improved operators have a better projection onto 
low mass states the corresponding correlation functions fall less steeply than those of 
the unblocked ones, hence the signal-to-noise ratio is improved, leading to considerably 
smaller statistical errors. 

In summary, we find that blocking has little or no effect in the Higgs phase, where the 
original local operators exhibit a rather good projection onto the ground state in each 
channel. In the confinement phase, on the contrary, blocking turns out to be necessary in 
order to obtain any useful signal at all. This is particularly pronounced in the 1 chan- 
nel. It was also demonstrated by using a large set of operators, that in the confinement 
phase for time separations up to ten timeslices one is typically still rather far away from 
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Figure 3: Three differently blocked W operators in the 1 channel, Higgs (left) and 
confinement (right) phase. 



the asymptotic region where all excited states have died out. This implies that the latter 
are rather light compared to the ones in the Higgs phase, as we shall see more explicitly 
in the next subsections. This is precisely what spoils an easy mass measurement in the 
confinement phase. We have increased blocking levels on each kind of operator until we 
could explicitly identify the ones with the best projection onto the lowest states. Thus 
we can be sure that we have found the optimal operators that can be constructed from 



T0[) by means of the blocking technique described in subsection 3.2 



4.2 Correlations of eigenstates 

Now we discuss the correlations of the eigenstates of the matrix correlators Cij(t) in- 
troduced in subsection p. 3| . Consider first the 1 channel. Employing a basis composed 
of the three operators used in Figure |3| we obtain, after diagonalisation, the three sets 
of effective W masses shown in Figure ^ for the Higgs and confinement phases. Com- 
paring with Figure 0, we conclude that this procedure has only slightly improved the 
projection onto the ground state. However, it has clearly separated off the excitations. 
Even though it is not always possible to identify extended plateaux for these excited 
states, one can nevertheless conclude from the comparison between Higgs and confine- 
ment phases that the gap between the ground state and the excitations is much larger 
in the former. Clearly, in order to obtain more quantitative information about the exci- 
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Figure 4: T/ie i/iree lowest eigenstates in the 1 channel, Higgs (left) and confinement 
(right) phases. 

tations, one would need to choose a larger basis of operators. We shall not pursue this 
possibility here, since in this channel our main interest is in the ground state. 

Next, consider the three lowest states in the ++ channel shown in Figure [|. In the 
Higgs region, the situation is rather simple with an isolated Higgs ground state and 
a large gap to excitations. Because the excited states are much higher in mass, their 
correlation functions fall rapidly, one loses the signal after a few timeslices, and it is 
difficult to identify well-defined excited states. Here one would also need to increase the 
basis of operators and to reduce the lattice spacing, in order to improve the situation. In 
contrast, in the confinement phase the diagonalisation has isolated three distinct states 
which were mixed previously. In Table Q the coefficients (cf. eq. (|2^)) with which the 
individual operators contribute to various eigenstates are shown. The labelling is such 
that $i denotes our best operator for the ground state, $2 the one for the first excited 
state, etc. 

According to this analysis the ground state in the confinement phase consists pre- 
dominantly of R- and L-contributions. The next state has contributions from all types 
of operators, with a dominance by R and L. As in the spin-one case, the gap between 
the lowest and first excited states is much smaller than in the Higgs phase. The separate 
plateau of the plaquette operators survives diagonalisation, representing a rather definite 
state. Table [I] shows that the plaquette operators indeed have practically no overlap with 
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Figure 5: TTie i/iree lowest eigenstates in the ++ channel, Higgs (left) and confinement 
(right) phases. 



the ground state $1. Conversely, the other ++ operators do not contribute to the state 
$3, which thus appears to be an object composed exclusively of gauge degrees of freedom 
and very little mixing with operators containing scalars. In the pure gauge theory this 
object would correspond to a glueball. It seems natural to interpret this state in the 
Higgs model as a ++ 'W-ball" , composed almost entirely of gauge bosons. As shown in 
Table [I], the basis of eigenstates contains two more states $6, $9 with almost exclusively 
plaquette contributions, thus appearing to be excitations of the state $3. Details of the 



spectrum of excited states will be presented in subsection EO 



4.3 Mass spectrum and finite-size analysis 

Now we proceed to presenting our complete set of results for the spectrum of the SU(2) 
Higgs model in three dimensions. We perform an analysis of finite-size effects and finally 
extrapolate our results to the continuum limit. 

In Tables and ^| we summarise the results on all lattices and for all values of (3a, Ph 
used in our calculation. All masses quoted in this section have been obtained by fitting 
the correlation functions to the functional form in eq. As has been demonstrated 
in subsection [O], our signals for the lowest states show quite pronounced plateaux. The 
situation is more difficult for the excited states. In Tables |^ and |3| we only record masses 
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$i = HI 


$ 2 = El 


$ 3 = H3 


$ 6 = H3* 


$ 9 = H3** 


01 


= P 3 


0.932(1) 


0.315(4) 


0.034(7) 


0.0007(54) 


0.0019(13) 


02 


= P 4 


0.9962(2) 


0.017(4) 


0.024(3) 


0.0020(53) 


0.0018(11) 


03 


= i? 5 


0.850(2) 


0.493(4) 


0.07(1) 


0.0029(50) 


0.0041(16) 


04 


= P 2 


0.071(3) 


0.11(1) 


0.615(3) 


0.64(1) 


0.435(8) 


05 


= P 3 


0.068(3) 


0.19(2) 


0.973(4) 


0.017(4) 


0.107(2) 


06 


= P 4 


0.036(3) 


0.13(1) 


0.631(4) 


0.594(9) 


0.476(8) 


07 


= L 44 


0.818(2) 


0.200(5) 


0.048(6) 


0.01(4) 


0.005(7) 


08 


= L 54 


0.661(3) 


0.562(4) 


0.05(1) 


0.004(28) 


0.012(5) 


09 


= L 65 


0.337(3) 


0.655(4) 


0.021(13) 


0.08(7) 


0.04(2) 



Table 1: Coefficients a^- as defined in eq. fllP ) of the operators used in the simulation for 
the three lowest ++ states in the confinement phase (j3c = 12, (3h = 0.3411, L 2 -T = AO 3 ). 
In the header, we also introduce the labelling for scalar states used below. 

for which we could identify a plateau of at least three timeslices extension in an effective 
mass plot. Those cases where the statistical errors of the correlation function were large, 
or where the overlap of the diagonalised operators onto the desired state was small, 
are marked by an asterisk. 



Pa 


Ph 


L 2 -T 


aM H i 


aM H3 


aM w 


12 


0.3418 


20 3 


0.1944(13)(1) 


1.270(16)(5)* 


0.624(5)(1) 






16 2 • 32 


0.1955(8)(9) 




0.625(2) 


9 


0.3450 


14 2 • 20 


0.2627(14)(4) 


1.90(5)* 


0.836(2)(1) 


7 


0.3488 


20 3 


0.348(2)(2) 




1.067(4)(3) 



Table 2: Mass estimates in the ++ and in the 1 channels in the Higgs phase. The first 
error is statistical, the second is an estimate of systematic effects. 

We investigate finite volume effects in detail for (3q = 12. Numerically, the infinite 
volume limit is reached when the change in a mass with increasing lattice size is smaller 
than the statistical errors. In order to avoid additional finite size studies for the smaller 
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Ph 


L 2 ■ 


T 


aM H1 


aM H2 


aM m 


aM w 


12 


0.3411 


40 3 




0.2903(15)(12) 


0.514(4)(1) 


0.563(5)(2) 


0.447(8) (3) 






32 3 




0.2885(25)(15) 


0.440(9) (4) 


0.527(13)(13) 


0.442(4) (5) 






26 3 




0.2813(22)(5) 


0.334(10)(4) 


0.544(12)(2) 


0.443(9)(2) 






20 3 




0.226(6) 


0.306(8)(3) 


0.509(11)(10) 


0.423(10)(4) 






16 2 


• 32 


0.1739(32)(19) 


0.247(12)(5) 


0.540(7)(3) 


0.422(8) (4) 






10 2 


■ 30 


0.121(4)(1) 


0.136(8)(1)* 




0.469(H)(1) 


9 


0.3438 


26 3 


0.387(2)(4) 


0.677(14)(3) 


0.772(10)(5) 


0.610(4)(2) 


7 


0.3467 


30 3 




0.510(2)(2) 


0.912(12)(4)* 


1.003(6)(2) 


0.799(8)(6) 






20 3 




0.512(3)(3) 


0.908(12)(18) 


0.997(20)(3) 


0.801(4)(2) 



Table 3: Mass estimates in the ++ and in the 1 channels in the symmetric phase. 
The first error is statistical, the second is an estimate of systematic effects. 

values of /3g, we take the required spatial length corresponding to the large volume limit 
of the lattice at (3g = 12 in units of the Higgs mass, MhiL, and scale it down to the 
lower /^-values. This way we ensure that the simulations at the smaller values of $q are 
done in the same physical volume as for (3g = 12. Strictly speaking, this procedure is 
only valid if the considered range of values for (3q is in the scaling region, an assumption 
which turns out to be satisfied rather well, as we shall see a posteriori. After infinite- 
volume masses have been determined for different /^-values they can be extrapolated to 
(3 G -> oo. 

The large-volume limit of the Higgs phase is rather easy to reach. Table |2] gives the 
Higgs and W boson masses in lattice units as measured on lattices with spatial lengths 
L = 16, L = 20 at (3g = 12. It is seen that for both states the masses on the two lattices 
are compatible within the statistical errors. 

Again the situation is much more difficult in the confinement phase, as is illustrated 
in Figure ||]. There are strong finite size effects for the lightest scalar state, which are 
only under control for lattices larger than L = 32 (1/L = 0.031). We estimate that the 
infinite-volume limit for the scalar ground state in the confinement phase is reached for 
aMn\L ~ 10. The vector boson mass, on the other hand, shows only little dependence 
on the volume. The IV-ball is just getting close to the large volume limit on a lattice with 
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L = 40 (1/L = 0.025). The most pronounced finite size effects of all states investigated 
are displayed by the intermediate state $2- We conclude that the ground state masses in 
both channels have reached the infinite volume limit, while for the excitations it would 
be desirable to go to larger lattices. In order to get an estimate for the situation in larger 
volumes it is instructive to consider (3g = 7, where a given lattice size (here we consider 
L = 30) corresponds to a larger volume in physical units than at {3q = 12. Since the 
mass of the lowest state $i is free of finite size effects at L = 40 the ratio of this mass 
at the two /?c-values may be used to scale lattice size and masses according to 

MM = 3^.7) ^:^:^ 52.^.12), 

.«,<& = 12,* = 52.9) = ^.T.L.SO)^^^). (27) 

The result of this scaling is shown as the open data points in Figure || Now the state $2 
also seems to have approached the large volume limit. However, since it is very close to 
the TV-ball at large volumes one expects some mixing between these states. Comparing 
our data for the coefficients from the lattices with L > 26 we find growing admixtures 
of plaquette operators to $2 with increasing volume, while the composition of the state 
$3 remains rather unchanged by the approaching $2- In order to be absolutely sure that 
$2 really represents an isolated state further investigations are required on larger lattices 
or at different parameter values, where $2 arid $3 might be more clearly separated. 

In ref. || it was stated that at (3g = 12 the results for the lowest ++ state were 
practically indistinguishable on lattices of size 30 3 and 20 3 . In contrast to this, we find a 
rather strong dependence of clMhi in this range of lattice sizes. In particular, our results 
on 32 3 and 20 3 are clearly incompatible. We ascribe this to a presumably incomplete 
isolation of the ground state in ref. || . In addition, we observe that on the 20 3 lattice our 
vector boson mass in the confinement phase is about 35%, and the scalar ground state 
about 25% lower than those reported in ||. 

We conclude that the construction of improved operators is an indispensable tool in 
the study of the mass spectrum of our model in the confinement phase. In view of this, it 
would be very interesting to apply this technique to mass calculations in the vicinity of 
the phase transition, and at higher values of A3/ p| corresponding to more realistic zero 
temperature Higgs masses. 

The finite volume effects that we have analysed so far are to do with the size of 
the spatial volume. There are, however, additional finite volume effects which have to 
do with the finite extent in time of the lattice. Of course the eigenvalues of the transfer 
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Figure 6: Finite volume study for the confinement phase at (3g = 12. Squares, circles and 
diamonds represent the three lowest ++ states, whereas triangles denote the lowest 1 



state. Open symbols indicate the data extrapolated from (3c = 7 according to eq. (|F/[j 
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matrix and (lattice) Hamiltonian, H, are not altered by varying T. However, what may 
change is the relationship between those eigenvalues and the exponents in the decays of 
our calculated correlation functions. For example, the fact that the values of our masses 
are with respect to the vacuum energy arises because our expectation values contain 
the partition function, Z = Tr{e~ HT }, as a normalisation factor and this factor will 
normally be dominated by the vacuum energy. If T is sufficiently small, however, then 
Z may receive significant contributions from excited states and the masses we calculate 
from our correlation functions may be shifted by the corresponding shift in the 'effective' 
vacuum energy. Exactly what the effect of this is going to be is a complicated matter 
since, on the one hand, similar contributions occur in the numerator of the correlation 
function and this may lead to a partial cancellation of this correction. On the other hand 
our scalar masses involve the subtraction of a vacuum expectation value of the operator, 
and this will also be affected. Nonetheless, although we cannot easily estimate where such 
effects may be important, we note that since the leading correction to Z is 0(e~ aMHlT ), 
we need to be concerned once aM^i is small. 

To obtain a quantitative control over this potential problem, we have taken our 10 2 
spatial lattice at @g = 12 (because it has the smallest value of aMm and we have repeated 
the calculations, with the same basis of operators (which, unusually, happened to be 6 in 
this case), for T = 20 and T = 12. We have extracted masses in the same way as on the 
T = 30 lattice and have found aM m = 0.125(3) for T = 20 and aM H1 = 0.111(3) for 
T = 12. Thus there are no finite-T effects within these small errors down to T = 20 and 
even at T = 12, where e~ aMHlT ~ 0.24, the shift in the extracted mass is only ~ 10%. 
At T = 20, aMu\T ~ 2.5 and this gives us a benchmark value for judging when we 
should be safe from such corrections. We find no significant effects, within errors, for our 
other states. Of course, these effects may be somewhat different in the different phases, 
and, to the extent that scaling is violated, at different j3c- However the volumes that we 
use for extracting our final masses have values of clMhiT so much larger than the above 
benchmark value that we saw no reason to repeat this analysis in those other cases. 

4.4 Higher excitations 

The diagonalisation procedure also enables us to compute masses of more highly excited 
states, which were not mentioned in Table § . These, however, are determined with much 
less accuracy, since the variational basis for these states is smaller. We nevertheless find 
it instructive to give a qualitative discussion of that part of the spectrum. Since the gap 
between the ground state and excited states in the Higgs phase is rather large for both 
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Pa 


Ph 


L 2 -T 


aM* w 


aM* H3 


aM* H % 


12 


0.3411 


40 3 


0.682(24) 


0.840(24) 


1.02(2) 






32 3 


0.622(19) 


0.804(10) 


0.983(16) 






26 3 


0.636(23) 


0.773(12) 


0.974(30) 


9 


0.3438 


26 3 


0.854(34) 






7 


0.3467 


30 3 


1.193(17) 










20 3 


1.219(21) 







Table 4: Mass estimates for the excitations of the vector boson and the W-ball in the 
symmetric phase. 

the Higgs and the vector boson (as can be seen, e.g. in Figures^ and [5]), we restrict our 
discussion to the symmetric phase. 

In Table|4| we present the results for those states where we felt confident enough to 
quote a mass estimate. For the iy-ball, the correlation functions of the first and second 
excited states were those where plaquette contributions were clearly dominant (see e.g. 
Table p] for the overlaps of operators $6 and $9 at (3q = 12, L 2 ■ T = 40 3 ). 

It is instructive to compare the mass estimates for the IV-ball and its excitations 
with the glueball spectrum in the pure gauge theory. In Table |5| the masses in lattice 
units of these states are compared with those of the corresponding ++ glueball and its 
first two excitations at f3c = 12 in three-dimensional pure SU(2) gauge theory f2TS| ]. 
The striking agreement between these states in the two theories indicates a remarkably 
complete decoupling of the pure gauge sector from the Higgs part in the SU(2) Higgs 
model. 





aM aM* aM** 


SU(2) pure gauge 


0.563(5) 0.805(8) 0.982(14) 


SU(2) Higgs 


0.563(5) 0.840(24) 1.02(2) 



Table 5: Comparison of ++ glueball and W-ball and their first two excitations at (3q = 
12. 
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The existence of a separate FT-ball which does not mix at all with other ++ states 
is rather unexpected in view of the coupling between scalar and link variables in the 
tree-level action, and this suppression of mixing must be of dynamical origin. It would 
be interesting to see whether this isolation of the pure gauge sector persists also in the 
1 channel. 



4.5 How certain is the ground state? 

Measuring and diagonalising the correlation matrices provided us with valuable insight 
into the excitation spectrum of the theory. What can we say about the existence of very 
light states? In the ++ channel we have a nearly complete projection onto the lowest 
state, and in the 1 case the projection onto the lowest state looks quite acceptable as 
well. In the last section it was demonstrated that our operator basis includes the optimal 
operators which can be obtained from the operator types QUID by means of the blocking 
techniques (IT5J) , QT5D . Although we do find the lowest 1 state to be about 30-40% 
smaller than that quoted in ||, our lowest masses are still much larger than the Higgs 
and vector boson masses predicted by the gap equations f|. 

Of course, we cannot strictly rule out the existence of a lighter state which may only 
show up at distances larger than those up to which we have a good signal. If there were 
such states, however, they would have to have a rather poor overlap with our operator 
basis. This can be made more quantitative as follows. All effective masses presented 
so far were obtained from (p6|) under the assumption that the corresponding correlation 
function is dominated by a single lowest state. Let us now assume that there is one lighter 
state in each channel such that our measured correlation functions would correspond to 
a superposition of two states, 

Ci(t) ~ Ai (e" aMlt + e~ aKh (T -* } ) + A) (e~ amt + e^" 1 ^) , (28) 

where am corresponds to the supposed very light mass and aM ; is of the size of the mass 
we extracted assuming a single exponential correlation function as in eq. (|25|) . Fixing the 
assumed light mass am to values motivated by the study in @], we try to fit our data for 
the low states by the correlation function (|28|). Some results are presented in Table |6|. 

The amplitude A is consistent with zero in all cases. Adding two standard deviations 
to A Q we get a bound at 90% CL for the ratios of the amplitudes. The square root of 
this ratio, which is given in the last column of the table, may serve as an estimate of the 
maximal matrix element that a lower state has with the corresponding eigenstate of our 
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am = 0.1 


aMi 




A 


y/(Ao + 2a) /Ai < 


$1 


0.287(4) 


1.810(18) 


-0.016(17) 


0.10 




0.515(8) 


0.341(2) 


0.0004(17) 


0.11 


$ 3 


0.567(8) 


1.031(8) 


0.003(5) 


0.11 


<S>wi 


0.449(14) 


0.881(52) 


0.0004(25) 


0.078 


am = 0.07 


aMi 


M 


A 


y/{Ao + 2a) /A { < 


$i 


0.288(3) 


1.805(12) 


-0.0095(97) 


0.074 




0.515(7) 


0.3411(3) 


0.0003(12) 


0.089 


$ 3 


0.566(7) 


1.031(8) 


0.0020(34) 


0.092 


<$>Wl 


0.448(12) 


0.880(50) 


0.0002(15) 


0.06 



Table 6: Results of double exponential fits with a fixed assumed light mass to the lowest 
states in the confinement phase (fie = 12, (3h = 0.3411, L 2 ■ T = AO 3 ). 

basis. This suggests that it is rather unlikely that significantly lighter states have been 
missed. 

A potential source of systematic errors in the reported values of ground state masses 
is the residual contamination of the correlation function by higher excitations. Of course, 
the blocking procedure in conjunction with our variational technique is designed to opti- 
mise the projection onto the ground state. In the vector channel, however, the plateaux 
set in at larger values of t, thus showing that the ground state does not dominate the 
correlation function at very early timeslices. 

In order to quantify this systematic error, we performed a double exponential fit 
similarly to eq. (|2q) . Here, however, the mass am was fixed to the mass estimate for the 
first excited state in either the scalar or the vector channel. 

Extending the fitting interval to earlier timeslices, we found that the double expo- 
nential fit does not change at all the mass of the lowest ++ state, thus confirming that 
a nearly perfect projection has been achieved. In the vector channel, the double expo- 
nential fit gave slightly lower results for the mass of the ground state. We found that 
the mass decreased by about 5%, but that none of the mass estimates using a double 
exponential fit were incompatible with the result using a single exponential. 
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We conclude that practically all contamination from higher states has been eli- 
minated in the scalar channel, while higher excitations might lead to an uncertainty 
of about 5% in the mass of the vector channel. 

4.6 The continuum limit 

Our next task is to extrapolate the lattice spectrum to the continuum by taking the 
limit (3g — > oo. The continuum limit is performed only for the lowest states where our 
results are accurate enough. We have taken our results at all three /^-values at the 
largest respective lattice sizes, which, as Figure^ shows, have reached, or are close to, 
the infinite volume limit. 




Figure 7: Continuum limit in the Higgs (left) and confinement (right) regions. Squares, 
circles and diamonds represent the three lowest ++ states, whereas triangles denote the 
lowest 1 state. Open symbols indicate the data extrapolated to 1/(3q = 0. 

In the confinement phase, the dimensionless combinations aM /3g/4 = Mjg\ were 
extrapolated linearly in 1/ (3q for the three lowest ++ and the lowest 1 states. In the 
Higgs phase, only the lowest scalar and vector states were extrapolated, since higher 
excited states could not be clearly identified at all three /^-values. in Figure |7| the data 
at the three /^-values are shown together with the extrapolated results. In addition, 
we extrapolated the dimensionless ratio aMn\j aMyj linearly in in both phases. 

Table [7| shows a summary of the continuum values of the individual masses and the mass 
ratio for both phases. 
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\ 3 /g2 = 0.0239 




Mhi/qI 


M m /gl 


M m /gi 


M w /gl 


M H1 /M W 


Higgs, 

fil/gl = -0.020 


mass 
X 2 /dof 


0.547(12) 
0.80 






1.91(3) 
1.25 


0.282(8) 
1.70 


Confinement, 
fi 2 3 /gl = 0.089 


mass 
X 2 /dof 


0.839(15) 
0.74 


1.47(4) 


1.60(4) 
0.42 


1.27(6) 
0.06 


0.655(30) 
0.58 



Table 7: Continuum values of the three lowest scalar and the lowest vector states, as 
well as the ratio Mhi/Mw in the Higgs and confinement phases. Since the extrapolation 
of Mh2/ 1 g\ was performed using only the data at (5q = 12, 7, we cannot quote x 2 /dof ; 
the error is a subjective estimate. 



5 Summary and Conclusions 



We have presented results for the mass spectrum of the continuum SU(2) Higgs model 
in three dimensions at selected points in the symmetric and broken phases of the model. 
In order to get reliable mass estimates, the use of improved lattice operators turned out 
to be crucial. This is of particular importance for the investigation of the possibility of 
very low-lying states of the kind predicted by the analytic approach presented in |4[]. 

Using our particular blocking procedure, we were able to increase the projection onto 
the ground state dramatically. In most cases in the scalar channel, we achieved projections 
of essentially 100%, whereas in the vector channel values for the overlap ranged between 
75-95%. Undoubtedly, with a more refined smearing or blocking procedure, one could 
improve the signal for the ground state in the 1 channel even further. We wish to 
emphasise the importance of a high projection onto the desired state, since otherwise 
the possible misidentification of plateaux in the effectice masses is a source of large 
systematic errors which are difficult to quantify. Due to our use of the blocking procedure, 
we observe quantitative differences in the masses of the lightest scalar and vector states 
on specific lattices in the symmetric phase compared to ref. |J. Furthermore, we observe 
strong finite-size effects in the ground state of the ++ channel in the symmetric phase. 

Within the framework of our calculation we find no evidence for very small masses 
in the scalar and vector channels in the symmetric phase, as predicted by 0. We wish to 
point out, however, that we considered correlations of gauge-invariant composite opera- 
tors, whereas the correlators of elementary fields used in the analytic approach in j|] are 
gauge- dependent. There are indications from the numerical work reported in [p7] that 
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correlations of these gauge-dependent operators indeed exhibit a signal corresponding to 
a very low effective mass of the gauge boson. This needs to be better understood. 

Our computation of masses of excited states confirms the existence of a dense spec- 
trum of states in the confinement phase of the model. This appears to be consistent with 
the picture that bound states constitute the particle content of the symmetric phase. A 
surprising result of our calculation is the existence of states that are composed almost 
entirely out of gauge degrees of freedom. This 'W-ball" and its excitations are almost 
identical in mass to their gluonic counterparts in the pure SU(2) gauge theory. We are 
thus led to conclude that the pure gauge sector in the SU(2) Higgs model approximately 
decouples from the scalar degrees of freedom, a phenomenon which must be of dynamical 
origin. 

We have shown in this paper that by using various refined calculational tools in 
lattice simulations of the SU(2) Higgs model, detailed information of the mass spectrum 
in the symmetric phase can be gained. This is important for the development of effective 
theories of the symmetric phase, which will serve to analyse the nature of the phase 
transition at very large Higgs masses [[U], |28|, and to describe the thermodynamics of 
the electroweak plasma in the high temperature symmetric phase in the early universe. 
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